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We investigate the systematics of fission barriers in superfieavy elements in tlie range Z = 108- 
120 and = 166-182. Results from two self-consistent models for nuclear structure, the relativistic 
mean-field (RMF) model as well as the non-relativistic Skyrme-Hartree-Fock approach are compared 
and discussed. We restrict ourselves to axially symmetric shapes, which provides an upper bound 
on static fission barriers. We benchmark the predictive power of the models examining the barriers 
and fission isomers of selected heavy actinide nuclei for which data are available. For both actinides 
and superheavy nuclei, the RMF model systematically predicts lower barriers than most Skyrme 
interactions. In particular the fission isomers are predicted too low by the RMF, which casts some 
doubt on recent predictions about superdeformed ground states of some superheavy nuclei. For 
the superheavy nuclei under investigation, fission barriers drop to small values around Z = 110, 
N = 180 and increase again for heavier systems. For most of the forces, there is no fission isomer 
for superheavy nuclei, as superdeformed states are in most cases found to be unstable with respect 
to octupole distortions. 

PACS numbers: 21.30.Fe, 21.60.Jz, 24.10.Jv, 27.90. -|-b 



I. INTRODUCTION 

The search for superheavy elements (SHE) has made 
exciting progress in the last few years P, Q, H- New, 
often more neutron-rich, isotopes of elements Z = 108 
0, Z = 110 '55, and Z = 112 have been reported as 
well as the synthesis of the new elements Z = 114 and 
Z = 116 0. At the same time earlier experiments were 
confirmed and more data for already existing isotopes 
collected . Together with the a-decay products of the 
newly synthesized nuclei the known region of superheavy 
elements has grown substantially. 

SHE are by definition those nuclei at the upper end of 
the chart of nuclei where quantum-mechanical shell ef- 
fects reverse the trend of decreasing - and, for the heav- 
ier ones, practically vanishing - liquid-drop fission bar- 
riers to produce significant stabilization. The lifetimes 
of the recently found SHE are many orders of magni- 
tude smaller than the early optimistic estimates [lOl [ll| . 
Additionally the systematics of fusion cross sections sug- 
gests that the extension of the chart of nuclides might be 
limited by the production mechanism of SHE, not their 
decay Q . With recent experiments heading for unknown 
territory, theory has to provide reliable predictions for 
the stability and the most accessible regions in the land- 
scape of nuclides. A crucial feature is here spontaneous 
fission which is characterized by the fission barrier. 

Not too much is known experimentally on the fission 
barriers of transfcrmium nuclei. Although their height is 
not known, barriers for ■^^^No and ■^^"'No are high enough 
to stabilize these isotopes against fission up to angular 
momentum 20 fl^ . An analysis of all available data for 



fusion and fission of ^112, ^^^\IA and ^^116 was given 
recently in Ref . (l^ . Surprisingly the barrier heights de- 
duced for these very heavy nuclei are similar, or even 
slightly larger, than the ones of actinide nuclei in the 
240p-|^ region. 

For some other superheavy nuclei it is known that the 
barrier is relatively small from the simple fact that fission 
is their preferred decay channel. This leads to another 
aspect of the stability against fission: the experimental 
identification of new superheavy nuclides is much simpler 
when a decay is the dominating decay channel. All recent 
new decay chains from Dubna end in a region of fissioning 
nuclei which prevents an extension of the known region of 
SHE to the "southeast" with current experimental tech- 
niques. 

The calculation of fission half-lives is a very demanding 
task, which became clear quite early |14|. First explo- 
rations of the potential energy surfaces of transfermium 
nuclei demonstrated already that triaxial and reflection- 
asymmetric degrees of freedom often greatly reduce the 
fission barrier. The collective mass provides the met- 
ric for the dynamical calculation of the fission half-lives. 
There is no published work so far that considers all ingre- 
dients using self-consistent models. First ambitious steps 
in that direction were taken for the calculation of the de- 
cay of the fission isomer into the ground state . The 
published fission half-lives from macroscopic-microscopic 
mean-field models also simplify the task considerably by 
restricting shape-degrees of freedom to axially symmetric 
ones and using a phenomenological parameterization of 
inertia parameters. The detailed potential energy land- 
scape is also an ingredient for estimates on the fusion 
cross section, although somewhat different paths have to 
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be considered asymptotically, see e.g. Refs. and fur- 
ther references therein. 

It is the aim of this paper to analyze the extrapo- 
lation of self-consistent mean-field models, namely the 
Skyrme-Hartree-Fock (SHF) approach and the relativis- 
tic mean-field (RMF) model, concerning large-amplitude 
deformation properties of SHE. To that end, we present a 
systematic survey of fission paths and barriers for a broad 
range of SHE, scanning the a-decay chains of even-even 
nuclei that are accessible with the current experimental 
techniques. The nuclei considered in this study cover the 
proton numbers 108 < Z < 120 and the corresponding 
neutron numbers 166 < N < 182. We first benchmark 
our models by calculating axial fission barriers and iso- 
meric states for heavy actinide nuclei for which experi- 
mental data are available, namely isotopes with proton 
numbers ranging from Z = 90 to 98. 

Section ^ summarizes the most important results of 
earlier studies relevant for our calculations. Section IIIII 
explains the theoretical and numerical methods which are 
used in our study. In section Hvl the results for actinides 
are discussed. Section |V] presents detailed results for de- 
formation energy curves in all considered superheavy nu- 
clei and for all the different forces. Section IVTI discusses 
the result in terms of key quantities such as deformation 
energies and barriers. Finally, section IVlII attempts to 
identify the underlying reasons for the different predic- 
tions among the forces and models. 



II. EARLIER CALCULATIONS 

Estimates of the stability of SHE against fission have 
a long history. First explorations of the potential land- 
scapes were made with phenomenological corrections 
for shell effects 0, ^J, even before Strutinsky intro- 
duced the microscopic-macroscopic method [18], which 
was then immediately applied to large-scale calculations 
of fission barriers of heavy and superheavy nuclei, see e.g. 
Ref. [l9l |. In this framework, the deformation energy is 
minimized in a limited space of shape parameters. It soon 
became clear that reflection-asymmetric ,2Qj and triaxial 
PH shape degrees of freedom have to be included. Com- 
plementing the collective potential energy surfaces by 
mass parameters enabled the calculation of fission half- 
lives of heavy and superheavy nuclei, with phenomeno- 
logical mass parameters |0, 12^ as well as with mi- 
croscopically computed cranking masses Studies 
along that line are continued until today with improved 
parameterizations of the mic-mac model, either the 
finite-range droplet (FRDM) plus folded- Yukawa single- 
particle potential model and phenomenological masses 
[24j or the Yukawa-plus-exponential (YPE) macroscopic 
plus a Woods-Saxon microscopic model and cranking 
masses j2^ |2^ Nearly all large-scale calculations 

of fission lifetimes, however, consider axially symmetric 
shapes only (as we will do). A recent exception is pre- 
sented in Ref. 12^ where triaxial shapes are taken into 



account, however, at the price of a reduced number of 
shape degrees of freedom in other places. 

There do exist also systematic calculations of fission 
barriers where the macroscopic part of the energy is cal- 
culated within the semi-classical Thomas-Fermi approx- 
imation [23. One step further toward self-consistency is 
the extended Thomas-Fermi Strutinsky-integral (ETFSI) 
approach where the same microscopic Skyrmc force is 
used to calculate the macroscopic part of the binding 
energy and to determine the single-particle spectra for 
the calculation of the shell correction. For a large-scale 
survey of the axially and reflection-symmetric potential 
energy surfaces of heavy and superheavy nuclei in the 
ETFSI approach using the Skyrme interaction SkSC4 see 
Refs. |2^|33|. Results allowing also for triaxial degrees 
of freedom are presented for selected nuclei in Ref. |31l| . 

To determine the fission path in microscopic- 
macroscopic, Thomas-Fermi, and ETFSI calculations, 
the total energy is minimized with respect to parame- 
ters of the nuclear shape. There are numerous parame- 
terizations to be found in the literature, which differ by 
emphasizing either high-order multipoles at small defor- 
mation, or the fragment deformation of two-center-type 
configurations at large deformations, see Ref. ,33] for an 
overview. While the latter are important for the proper 
description of the saddle-point of actinide nuclei (which 
is located at very large P2), the first might be better 
suited for superheavy nuclei where the saddle point is lo- 
cated at rather small ^2- The number of shape degrees 
of freedom in actual calculations does rarely exceed five, 
in most cases it is even smaller. 

This restriction does not exist in the framework of self- 
consistent models. Besides very general spatial symme- 
tries that are imposed (e.g. axiality, reflection symme- 
try, or triaxiality) there are no further assumptions made 
on the nuclear density distribution. As this makes self- 
consistent calculations more costly in terms of computa- 
tional time, there is much less published work employing 
self-consistent models so far. 

A systematic study of the deformation energy of super- 
heavy nuclei along the valley of /3 stability in the region 
100 < Z < 128 and 150 < TV < 218 in HFB calculations 
with the Gogny force Dls under restriction to axially 
and reflection symmetric shapes was presented in Ref. 
[33I ]. Potential energy surfaces of selected heavier nu- 
clei are presented in Ref. [s^]. The full potential energy 
surface in the /?-7 plane of a few selected nuclei as re- 
sulting from SHF calculations in a triaxial representation 
is presented in Ref. [35^ . This investigation points out 
the importance of triaxial shapes at small deformations 
(32 < 0.6. The fission barrier of some superheavy nuclei 
is reduced to half its value when relaxing the constraint 
on axial symmetry and going through triaxial paths. The 
fission path returns to axially symmetric shapes at larger 
deformations. But here it is necessary to allow for re- 
flection asymmetric shapes to accommodate the usually 
asymmetric fission. A first exploration of asymmetric 
shapes of SHE within self-consistent SHF and RMF mod- 
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els was presented in Ref. [Sg- For many superheavy 
systems, reflection-asymmetric shapes lower the fission 
path at large prolate deformation /32 > 0.6, and remove in 
most cases the outer barrier known from actinide nuclei 
(and persisting for superheavy nuclei when considering 
reflection-symmetric shapes only). 



III. THEORETICAL FRAMEWORK 
A. Effective Interactions 

We explore the potential landscapes using two widely 
used self-consistent mean-field models, namely the non- 
relativistic Skyrme-Hartree-Fock (SHF) method as well 
as the relativistic mean- field (RMF) approach 37]. There 
exists a great variety of parameterizations for both mod- 
els which often differ when extrapolated. It is long known 
that different parameterizations of a self-consistent model 
are not equivalent for the calculation of fission barriers 
of actinide nuclei, see Re f. l39| for a comparison of 
early Skyrme forces and |43| for a comparison of RMF 
forces. Comparisons between SHF and RMF hint at gen- 
uine model differences "i^, "ill. Relevant for our study 
is also that there exist conflicting predictions for the lo- 
cation of the spherical magic numbers in SHE, see Ref. 

E3i which can be expected to be reflected in the 
structure of the potential landscapes. 

A fair survey of the extrapolation of the models to 
large mass number and large deformation has therefore 
to cover a selection of typical parameterizations. We have 
chosen parameterizations which give a very satisfactory 
description of stable nuclei but differ in details. Namely 
we use the Skyrme interactions SkP 0, SLy6 [S], SkI3, 
and SkI4 In the relativistic calculations the param- 
eterizations NL3 113 and NL-Z2 43] of the standard La- 
grangian are employed. 

The parameterization SkP has the isoscalar effective 
mass mg/m = 1 and was originally designed to describe 
the particle-hole and particle-particle channel of the ef- 
fective interaction simultaneously (we do not make use of 
this particular feature). The forces SLy6, SkI3, and SkI4 
stem from recent fits including already data on exotic 
nuclei (all three forces) and even neutron matter (SLy6). 
Both SkP and SLy6 use the standard spin-orbit interac- 
tion. The forces SkI3/4 employ a spin-orbit force with 
modified isovector dependence. SkI3 contains a fixed 
isovector part analogous to the non-relativistic limit of 
the RMF, whereas SkI4 is adjusted allowing free vari- 
ation of the isovector spin-orbit term. The RMF force 
NL-Z2 is fitted in the same way as SkI3 and SkI4 to a 
similar set of observables. 

A quantity that characterize the average deformation 
properties of an effective interaction is the surface en- 
ergy coefficient. It is determined for the model system 
of semi-infinite nuclear matter, which offers the cleanest 
procedure to define a surface energy, see, e.g. |33,EEE3 
and the references given therein. In table we give 



TABLE I: Compilation of bulk properties for the parame- 
terizations employed in tiiis study. The upper block shows 
the volume parameters incompressibility modulus K, effec- 
tive mass rriQ/m, and asymmetry energy coefficient ttaym. The 
lower block shows the surface energy coefficient Ogurf obtained 
from semi-infinite nuclear matter calculations. For the RMF, 
where the effective mass is momentum-dependent, mj/m is 
given at the Fermi momentum k = kp. This value, which is 
about 10% larger than the often quoted value for k = 0, de- 
termines the average level density around the Fermi energy, 
see also Ref. jij and references therein. 



force 


SkP SkI3 SkI4 SLy6 NL-Z2 NL3 


K (MeV) 


202 258 248 230 


172 


270 


mo/m 


1.00 0.59 0.65 0.69 


0.64 


0.67 


a^yra (MeV) 30.0 34.8 29.5 32.0 


39.0 


37.4 


asurf (MeV) 18.2 18.3 18.3 17.7 


17.7 


18.5 



values for bulk and surface properties of nuclear matter 
obtained exclusively in fully self-consistent Hartree-Fock 
calculations of semi-infinite matter js^, E| • Often values 
obtained within the extended Thomas-Fermi approxima- 
tion are given |39l ]. which are usually smaller by about 1 
MeV. 

A correction of the binding energy for spurious 
center-of-mass motion is performed as usual. For 
SkI3, SkI4, SLy6, and NL-Z2 the cm. correction 
Ec.m. = (Pc.m.)/2wA is subtracted after variation. For 
SkP the diagonal part of i?c.m. only is considered before 
variation, while for NL3 the harmonic-oscillator estimate 
for i?c.m. is subtracted, see Ref. details. The cm. 

correction, howe ver^ varies only little with deformation, 
see e.g. Ref. 0,153 • Note that the various recipes for 
cm. correction cannot be easily interchanged as their dif- 
ferences are partially absorbed into the force parameters. 
This has a consequence relevant for our study: The differ- 
ence between the "exact" and the approximate schemes 
used for SkP and NL3 scales in leading order as ~ A^^^ . 
During the fit this difference is incorporated into the ef- 
fective interaction which leads to the significantly larger 
surface tension for SkP and NL3 found in tabled see Ref. 
|l4h] . The deformation energy from otherwise equally fit- 
ted forces might differ on the order of 5 MeV at the outer 
barrier in actinides. For SHE, where the saddle-point is 
at smaller deformation, this effect can be expected to be 
less pronounced, but still might cause differences of a few 
MeV between forces. 

We treat pairing correlations within the BCS approxi- 
mation using an effective density-independent zero-range 
delta pairing force with the strength adjusted for each 
mean-field parameterization separately as described in 
Ref. .53]. Including a density dependence of the effective 
pairing interaction or an approximate particle-number 
projection might alter the barrier heights on the order 
of 1 MeV for actinide nuclei 'Hi'] . 

The coupled mean-field equations for both SHF and 
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RMF models are represented on a grid in coordinate 
space using a Fourier representation of the derivatives 
and are solved with the damped gradient iteration 
method as described in Ref. [s^. The numerical codes 
for both models |54, _56] share the same basic numerical 
routines which allows for a direct comparison of the re- 
sults. Note that the accuracy of grid techniques is fairly 
independent on deformation, which is an advantage to 
calculations using an harmonic oscillator basis expansion, 
see e.g. Ref. [S^l for a convergence study for ^^^Fm. 

Finally, beyond-mean-field effects can modify the fis- 
sion barrier. The most important corrections to the bind- 
ing energy remove the contributions from spurious collec- 
tive vibrational and rotational states which are inevitably 
admixed to the mean-field wave functions IH^ IfiTf . 
These corrections lower the barriers, typically up to 1 
MeV for the inner one and up to 2 MeV for the outer 
when starting from a well-deformed prolate ground state. 
The situation is less clear for transitional or spherical nu- 
clei. 

Altogether, the present study has uncertainties on the 
deformation energy for given configuration of the order 
1-2 MeV. Most of the possible improvements increase the 
binding energy. Thus we can assume to explore an upper 
limit for the barriers. The comparison of barriers between 
different forces is more robust because most corrections 
can be expected to be similar for all forces. 



B. Shape Degrees of Freedom 

We will consider reflection-symmetric as well as 
reflection-asymmetric fission paths. But we restrict the 
considerations to axially symmetric shapes. Our investi- 
gation also covers the prolate fission path only. It has to 
be kept in mind that novel fission paths may emerge for 
the heaviest of the nuclides discussed here, which start 
out from strongly oblate shapes and proceed through tri- 
axial deformations 2Lj. For this and the reasons given 
above our results provide an upper limit for the (static) 
fission barriers. This limitation holds also for most other 
work using self-consistent models published so far, as well 
as most of the results from mic-mac approaches. 

The deformation energy curves are obtained with a 
constraint on the mass quadrupole moment Q20 — {Q2o)- 
For reflection-asymmetric shapes, we also fix the center- 
of-mass with a constraint on the mass dipole moment 
(Qio) = 0. The constraints are added to the energy func- 
tional by means of Lagrange multipliers Besides 
these constraints, the deformation energy is minimized 
with respect to all axial multipole moments Qio for pro- 
tons and neutrons separately. In a self-consistent calcu- 
lation, the energy is not only minimized with respect to 
the deformation, but also the radial profile of the den- 
sity distribution, again separately for protons and neu- 
trons. With that, a self-consistent calculation explores 
many more degrees of freedom than the best microscopic- 
macroscopic calculation available so far. Some conse- 



quences will be discussed in section fVI CI below. 

The deformation energy is shown versus the dimension- 
less multipole deformations of the mass density which are 
defined as 



Pi 



in 



(/ Yio) with To = 1.2 fm. (1) 



Note that the (3e are computed from the expectation val- 
ues of the actual shapes and need to be distinguished 
from the generating moments which are used in mul- 
tipole expansions of the nuclear shape in microscopic- 
macroscopic models fs^. 

The constrained calculation does not always follow ex- 
actly the static fission path, which is defined as the path 
that follows the steepest descent in the multidimensional 
energy surface. Instead, for each value of Q20, one ob- 
tains a state which corresponds to a minimum with re- 
spect to all other degrees of freedom |5^. This might 
cause some problems to keep track of the path whenever 
the fission path has a small component only in the di- 
rection of the constraint. Often there exist two or even 
more distinct valleys in the multidimensional potential 
landscape which are separated by potential barriers. De- 
pending on the choice for the initial wave functions, the 
constrained calculation might find the nearest relative 
minimum only, which is not necessarily the absolute min- 
imum for a given constraint. The existence of distinct 
valleys complicates the interpretation of the deformation 
energy curves. When the solution jumps from one val- 
ley to another it misses the saddle point in between. The 
resulting uncertainty is not clear a priori, as only a calcu- 
lation including two or even more constraints can clarify 
if there is a flat plateau or a mountain ridge between 
the two valleys. The change from one to another valley 
in the potential landscape is accompanied by discontinu- 
ities in higher multipole deformations which can be used 
to identify them. In some cases the existence of two dis- 
tinct valleys might be the artifact of imposed symmetries, 
c.f. the case of ^^^Fm, where the two distinct valleys ob- 
tained from axially symmetric calculations are spurious 
as the two solutions are smoothly connected through tri- 
axial shapes [3^. 



IV. BARRIERS IN ACTINIDE NUCLEI 

Actinides are the heaviest nuclear systems for which 
data on the structure of the fission barrier are available. 
We use these nuclei to benchmark our models and forces, 
and to examine the force dependence of the predictions. 
We confine our investigation to axial barriers but release 
reflection symmetry in the calculations of the outer bar- 
riers and isomeric states. 

There exists a wealth of information about the (in most 
cases) double-humped fission barriers of actinide nuclei, 
see Ref. |59j and references given therein. The generic 
features of the static fission path are shown in figure ^ 
for the example of ^''"Pu calculated with the Skyrme force 
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FIG. 1: Example for the double- humped fission barrier of the 
typical actinide nucleus ^'"^Pu. The dotted line denotes an ax- 
ial and reflection-symmetric calculation, the full line denotes 
a triaxial (inner barrier) and axial and reflection-asymmetric 
calculation (outer barrier). The various shapes along the ax- 
ial paths are indicated by the contours of the total density at 
po = 0.07 fm-^ 



SLy6. The deformed ground-state has a calculated defor- 
mation of P2 = 0.29, which is in perfect agreement with 
the value of j32 — 0.29 that can be deduced within the 
rigid rotor model from the B{E2) t value of 13.33 ± 0.18 
obtained from Coulomb excitation [b^- The de- 
formation energy of 15.9 MeV of the ground state cor- 
responds to 0.9 % of the total binding energy. The in- 
ner barrier explores triaxial degrees of freedom, which 
reduce the (axial) barrier by about 3 MeV. There is a su- 
perdeformed fission isomer at f32 ~ 0.8 at an excitation 
energy of 3.0 MeV, which is somewhat larger than the 
experimental value of 2.25 ± 0.20 MeV Jlj. The outer 
barrier explores reflection-asymmetric shapes. The po- 
tential landscapes of adjacent actinide nuclei are similar, 
though for some nuclides there might appear a second 
isomeric state |6^ . 

The double-humped fission barrier of ^'^^Pu has often 
served as a benchmark for mean-field models, see Ref. 
[63l | for results obtained using Skyrme interactions, 
usin g Go gny forces and Refs. jifl, |63 for the RMF, see 
also |4l|. Early comparisons of barriers obtained with 
different Skyrme forces using approximations to full self- 
consistency were published in Refs. jsl, |3^ . A direct 
comparison of the potential energy curves with triaxial 
inner and reflection-asymmetric outer barriers obtained 
with the Skyrme interactions SLy6, SkM* and SkI4, the 
Gogny force Dls and the RMF forces NL3 and NL-Z2 
is presented in ,,37\ . The excitation energy of the flssion 
isomer has been studied with a variety of Skyrme forces 
in Refs. 

The influence of correlations on the deformation energy 
in the framework of Skyrme mean-field calculations has 
been studied in Ref. j69j. An exact angular momentum 
projection lowers the axial inner barrier by a bit less than 
1 MeV, the excitation energy of the isomeric minimum by 
about 1 MeV, and the reflection-symmetric outer barrier 



by about 1.5-2 MeV. Removing spurious quadrupole vi- 
brations from the mean-field states by configuration mix- 
ing of the angular-momentum projected mean field states 
lowers both the ground state and the isomeric state by a 
few 100 kcV. For ^'''"Pu, this effect is more pronounced 
for the isomer, which lowers its excitation energy even 
further, but it cannot be expected that this will be the 
same for all actinides. 

In this section we will confine ourselves to the heights 
of the inner and outer barrier as well as the excitation en- 
ergies of the isomeric states while postponing a thorough 
discussion of the potential landscapes to future work. As 
our goal is the extrapolation of the models to superheavy 
nuclei, a discussion of these key quantities suffices. 

The selection of even-even nuclei for this study is 

230rpv, 234rpi^ 234tt 238tt 238t3,, 242p-,, 

90 -Ltll40j 90 -Lni44, 33 U142, 92 Ui46, 94 -t^Ui44, 94 rUi48, 

if Pui52, if Cmi46, if Cmiso, if Cmi54, ifCfiss, which 
is only every second known nucleus in neutron number. 
We omit Fm isotopes as there are two competing paths 
at the outer barrier, and there is no continuous axial path 
for the inner barrier for some Skyrme forces |36j |. 

Figure|21compares calculated and experimental heights 
of the inner barrier. In the case of the inner barrier, the 
comparison of data with values obtained from axial calcu- 
lations is somewhat dangerous, as the static inner barrier 
is known to be triaxial and the energy gain through triax- 
ial deformation may differ for each force. Still, there are 
several conclusions that can be safely drawn from Fig.|21 
Our selection of forces suggests that there is a difference 
between SHF and RMF models. All Skyrme forces pre- 
dict that the inner barrier increases with neutron number 
up to iV = 150, most pronounced for SLy6. The Z depen- 
dence of the barrier height is most pronounced for SkI3 
and SkM, the Skyrme forces with an extended spin-orbit 
interaction, while it is negligible for SkP, the only force 
in our sample with the large effective mass mg/m = 1.0, 
c.f. table m One might suspect that the corresponding 
large level density suppresses shell effects compared to 
the other forces. On the other hand, for the two RMF 
forces the barriers stay nearly constant with A'', and show 
an increase with Z only. This finding suggests a signif- 
icant difference in shell structure between the SHF and 
RMF models. Experimental data do not show any sig- 
nificant dependence on A^ or Z at all, they just fall off a 
few 100 keV with mass number. From our present cal- 
culations, it cannot be decided if adding triaxiality will 
give a similar trend. 

There is also a difference in absolute height between 
SHF and RMF. With the exception of SkP, the inner 
barriers from SHF are significantly higher by about 2 
MeV than those obtained within the RMF. This is re- 
flected in the values for the mean deviation of the inner 
barrier height (1/?^) 1^-^*1 from the experimental 

value in MeV is 4.5 (SLy6), 2.7 (SkP), 4.1 (SkI3), 4.1 
(SkM), 1.5 (NL-Z2), and 2.3 (NL3) respectively, although 
this quantity is of limited signiflcance. Angular momen- 
tum projection will lower the barrier by about 1 MeV, 
and one can speculate only about the effect of triaxiality 
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FIG. 2: Height of the inner barrier from axial and reflection- 
symmetric calculations. Th {Z = 90), U (Z = 92), Pu 
(Z = 94), Cu (Z = 96), and Cf (Z = 98) isotopes are de- 
noted by open circles (for N = 140, 142), open triangles, open 
squares, open diamonds, and again open circles (A'^ = 152), re- 
spectively. Experimental data (full symbols) are taken from 
Data points for the same element are connected by lines. 



and other correlations on the trends with N and Z. It is 
tempting to assume that, when including these missing 
corrections, the SHF might still overestimate the barri- 
ers of the heavier nuclei but be on the right order for 
the lighter ones, while the RMF will underestimate the 
barriers of Th isotopes already on the mean-field level. 

The outer barrier heights as predicted by the various 
mean-field forces are shown in Fig. |3| For the outer bar- 
rier, it can be expected that our calculations cover all 
necessary degrees of freedom, so data and calculated val- 
ues can be directly compared. There are differences in 
absolute barrier height. With a mean deviation of the 
outer barrier height in MeV of 2.2 (SLy6), 4.7 (SkP), 1.4 
(SkI3), 3.8 (SkI4), 1.2 (NL-Z2), and 0.8 (NL3) respec- 
tively, SkI3, NL-Z2 and in particular NL3 give a quite 
good description of the barrier height on the mean-field 
level. The differences found for the barrier height seem 
to manifest themselves mainly in an overall offset, while 
all models and forces predict quite similar trends of the 
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FIG. 3: The same as Figure |51 but for the outer barrier. 



barrier heights with N and Z , a bit more pronounced for 
Skyrme forces, and a bit more damped for NL-Z2. The 
good news is that this overall trend is quite close to the 
experimental findings. 

As barriers from NL-Z2 are already always smaller 
than the experimental values, there is no room left for 
correlation effects. For shapes around the outer barrier, 
the (missing) rotational correction can be expected to 
be about 1.5 MeV larger than for the ground state of a 
well-deformed nucleus. When removing those 1.5 MeV 
from the outer barrier heights shown in Fig. 13 the bar- 
riers from the RMF are too low, in particular for NL-Z2, 
where not much will be left. 

The excitation energy of the fission isomer is displayed 
in Fig. ^ We have added also results obtained for the 
nuclei with neutron numbers between those used to in- 
vestigate the barriers above. To the best of our knowl- 
edge, only three experimental values for superdeformed 
O''^ states in even-even nuclei are available so far from 
spectroscopy in the superdeformed and normal-deformed 
wells, which are ^36,8^ g^j^^ 240p^^ ^^le recent collec- 
tion of data in • There are more superdeformed levels 
known in some other adjacent even and odd nuclei, but 
their quantum numbers could not be established so far. 
They all have in common that their excitation energy 
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FIG. 4: Excitation energy of the isomer, obtained from axial 
calculations. 



is at least 2 MeV, which also sets some constraints to 
our calculations. These values are also consistent with 
the data obtained from fits to fission-isomcr excitation 
functions, see, e.g., Ref. fflf. 

Although the sparse data do not allow for a detailed 
analysis, there are a few conclusions that can be drawn. 
On the mean-field level, SkP gives rather high energies 
and overestimates the data by about 2 MeV. Values from 
SLy6 are scattered around the data. RMF models pre- 
dict very low excitation energies. The Skyrme interac- 
tions SkI3 and SkI4 with extended spin-orbit interactions 
also underestimate the known excitation energies at least 
for certain elements. It is noteworthy that SkI3, NL3 
and particularly NL-Z2 predict the superdeformed state 
to be the ground state for some actinide nuclei, in con- 
tradiction with experimental knowledge about the spec- 
troscop y an d decay of those nuclei. It was already noticed 
in Ref. [i^l for selected examples that the RMF underes- 
timates the excitation energy of the fission isomer. This 
seems to be a general shortcoming of the RMF model, at 
least of most, if not all, of its standard parameterizations. 
This finding is not restricted to actinide nuclei, but was 
also observed for superdeformed states in the neutron- 
deficient A « 190 region [U Iz^, where it can be cured 



to some extend taking additional information about the 
spherical shell structure of ^°*Pb into account during the 
fit of the force parameters [t^ ■ The overbinding of the 
fission isomer can be expected to become even more pro- 
nounced when corrections for breaking of rotational and 
other symmetries are considered. 

From our small selection of Skyrme forces, it is hard 
to disentangle the infiuence of the effective mass and of 
the various spin-orbit functionals from the influence of 
the actual fitting procedure on the predictions for barrier 
heights. It is tempting to correlate the large deformation 
energy obtained with SkP with its large effective mass, 
but this is ruled out by recent Skyrme-HFB mass fits 
with effective mass around 1.0 that deliver much smaller 
barriers ^ than SkP. Since NL-Z2 and NL3 have the 
same functional form, the tendency of NL3 to larger bar- 
riers must result from the different strategies to adjust 
the force. 

It is a bit surprising that the differences in surface 
tension visible in table are not reflected in the inner 
barrier heights and cannot be solely responsible for the 
differences obtained for outer barrier heights and flssion 
isomers. Shell effects seem to play a much larger role for 
these quantities than the nuclear matter properties. For 
the forces in our sample, the difference in surface ten- 
sion seems to be compensated by other features of the 
forces through the flt. However, when the influence of 
shell effects is suppressed by comparing predictions of 
pairs of otherwise identically fitted forces with different 
surface tension, one finds indeed the expected difference, 
see Ref. the example of the Skyrme interactions 

SLy4 and SLy6 and the RMF forces NLl and NL-Z. 

Summarizing, a differences in the models, can be rec- 
ognized: SHF gives usually higher barriers than RMF. 
In comparison with experimental data, it seems that the 
RMF predictions are too low for barriers and fission iso- 
mers. These findings were already hinted in earlier inves- 
tigations of actinide [iol I4H and superheavy nuclei |3^ , 
but emerge even more clearly for the present systematic 
investigation of nuclei. 

In the next section we will see how these trends trans- 
late to superheavy nuclei. 

V. POTENTIAL LANDSCAPES 

Figure [S] provides a summary view of the deformation 
energy curves along the fission paths for all SHE under 
consideration here for the Skyrme interaction SLy6 and 
the relativistic mean-field force NL-Z2. The full lines de- 
note the asymmetric fission path, the dashed lines the 
symmetric fission path (which coincide at small defor- 
mations). The global trends are common for both forces 
(and also the others employed in this study) and in qual- 
itative agreement with earlier studies in mic-mac and 
semiclassical models: 

• There is a gradual transition from well-deformed 
nuclei with (32 ~ 0.3 around the deformed Z — 108 
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and N = 162 shell closures to spherical shapes ap- 
proaching N = 184. Note that earlier studies sug- 
gest that the neutron number is more important 
than the proton number to determine the ground- 
state shape. 



• Intermediate systems around Z = 114, N = 174 
have two distinct prolate and oblate minima at 
small deformation. On the basis of our axial cal- 
culations it is not clear if this leads to shape co- 
existence in these nuclei as concluded on similar 
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grounds in Refs. [7^|2Elz3- Calculations including 
triaxial degrees of freedom suggest that for some, 
perhaps all of these systems the prolate and oblate 
"minima" are connected through triaxial shapes 
without a barrier ,7S] . 

• While the ground-state deformation moves to 
smaller values, also the saddle point is shifted to 
smaller deformations. This means that the width 
of the fission barrier is not necessarily larger for 
nuclei with spherical shell closures. 

• The transitional nuclei in between belong to a 
regime of low fission barriers. This is consistent 
with the current (still sparse) experimental knowl- 
edge. The recent data from Dubna interpreted as 
the a-decay chain of ^^^116 indicate that fission is 
the preferred decay channel of ^®°110i7o, but not 
the heavier nuclides in this chain |^. 

• Above Z — 108, the static fission path switches 
from symmetric to asymmetric fission. 

• For most of the nuclei above Z > 110 considered 
here, reflection-asymmetric shape degrees of free- 
dom remove completely the outer barrier that is 
well-known from actinide nuclei, leading to a single- 
humped fission barrier only. 

• The superdeformed minima around /32 ~ 0.5 ob- 
tained from reflection-symmetric calculations with 
NL-Z2 for 298ll8i76 in Ref. |7| and ^921 16175, 
288ll4i74, '*^112i72, and 280ll0i72 in Ref. ^ are 
not stable with respect to octupole distortions, 
which makes the conclusions about superdeformed 
gro und states of superheavy nuclei drawn in Ref. 
|76l questionable. This finding is not com- 
pletely general as for some forces there remains a 
very small asymmetric barrier. 

The general features of these trends are easily understood 
in the more intuitive language of the mic-mac models (al- 
though they apply, of course, to the self-consistent mod- 
els as well). For nuclei at the lower end of the region 
investigated here, the potential energy surface from the 
LDM is rather flat around the spherical point and drops 
off fast at prolate deformations about ^2 ~ 0.7, c.f. fig- 
ure 7 in Ref. [s^. The structures seen in Figures 13 are 
mainly determined by the variation of the shell correction 
with deformation. The maximum of the shell correction 
follows the shell closures from the deformed Z = 108 and 
N = 162 to the spherical iV = 184 sheU. The potential 
wells are deepest in the vicinity of closed shells. With 
increasing Z, the plateau is shifted toward oblate defor- 
mations, while the LDM surface drops at smaller and 
smaller prolate deformations, which cannot be counter- 
weighted by the variation of -Eghcii- With that the saddle 
point moves in toward smaller deformations. For larger 
systems than those discussed here the potential energy 
surface becomes also unstable on the oblate side. 




/52 /32 

FIG. 6: The fission barriers of four selected nuclei normalized 
to the ground-state for four mean-field parameterizations as 
indicated. Shown is always the energetically favoured fission 
barrier. 



Figure Ogives also an idea where triaxiality might play 
a significant role: Whenever the deformation energy is 
smaller on the oblate side than for the same deforma- 
tion on the prolate side outside the prolate saddle point, 
oblate shapes might be unstable through a triaxial path. 
Of course this is neither a sufficient nor necessary condi- 
tion. 

While there is overall qualitative agreement among the 
two forces (and models), there are significant differences 
at a quantitative level. The RMF force NL-Z2 predicts 
lower barriers when going toward heavy systems than 
the Skyrme interaction SLy6. This is confirmed when 
directly comparing the deformation energy for selected 
nuclei, see Fig. El iggHsios is a well-deformed nucleus 
located at the edge of the "rock of stability" around 
fgjHsios, 182I2O is close to the spherical neutron shell 
TV = 184, while the other two are located in the transi- 
tional region. The figure compares more forces, now two 
from SHF (SLy6 and SkI3) and two from RMF (NL3 and 
NL-Z2). 

There are two kinds of differences: first, the systematic 
difference between SHF and RMF models which we saw 
already for the actinides, with the RMF giving smaller 
barriers, persists to superheavy systems, and second, an 
additional difference between the two RMF parameteri- 
zations NL-Z2 and NL3 concerning the outer barrier oc- 
curs, with NL3 being the only force predicting a double- 
humped barrier for the heavier systems. 
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FIG. 7: Deformation energy for the lowest minimum at small 
deformation (/32 < 0.4) for the nuclides and forces as indi- 
cated. White squares denote deformation energies smaller 
than 0.5 MeV. See the main text for a detailed description. 



FIG. 8: Height of the outer, usually asymmetric, fission bar- 
rier with respect to the isomeric state. White squares indicate 
a second bump smaller than 0.5 MeV. 



Comparing the potential energy curves from NL-Z2 
and SLy6 at small deformation for the heaviest nuclei, 
one sees also some differences concerning how strongly 
the nuclei are driven to sphericity. Comparing nuclides in 
the "northeastern" corner of figure [3 the onset of spher- 
ical ground states is predicted to be earlier with NL-Z2 
than with SLy6. This reflects the different predictions 
for shell closures from both models |4^ l4^ . While NL- 
Z2 (like all other standard RMF forces) predicts strong 
Z = 120 and N = 172 shells and a weak N = 184 sheU, 
SLy6 gives a strong TV = 184 shell and a weak Z = 120 
shell, which is not sufficient to guarantee a spherical 
ground state of this nucleus for non-magic neutron num- 
ber. Therefore nuclei at the upper end of figure [S] are 
much more driven to sphericity when calculated with NL- 
Z2 than with SLy6. 



VI. CHARACTERISTIC QUANTITIES 

The basic features of all the deformation energy sur- 
faces shown in the previous plots can be characterized by 
a few key numbers, i.e. the ground-state deformation en- 
ergy, the height of the inner fission barrier and the outer 
barrier, which will be discussed in this section. 



A. Ground-State Deformation Energy 

The ground-state deformation energy, i.e. the energy 
difference between the spherical shape and the (possibly 
deformed) ground state is plotted in figure [7| We con- 



sider only small deformations /?2 < 0.4. Zero deformation 
energy always indicates a spherical ground state, as we 
do not find coexisting well-deformed minima at the same 
energy as the spherical configuration. 

We see that all models and parameterizations pre- 
dict a transition from deformed nuclei around Z = 108, 
N — 162 to spherical N = 184 in agreement with earlier 
studies IzS • These deformation energies correspond 
to strongly prolate deformations at the lower corner of 
our selection of nuclei. For higher Z values, a shape iso- 
merism is established with two minima on the prolate 
and oblate side having approximately the same energy 
and only small deformations of (32 ~ ±0.15. Thus far the 
general trends agree. There are differences in quantita- 
tive detail, most prominently the fact that RMF predicts 
systematically lower deformation energies than SHF. 

It is to be noted that in SHE the neutron number most 
often determines the ground-state shape, while a magic 
proton number might not prevent deformation. This 
seems to be a general feature of all self-consistent models 
|35l l79| . and can also be observed in mic-mac models, 
which predict deformed Z — Hi isotopes 113,1^ far off 

= 184. This is due partially to the overall larger shell 
correction energy of the spherical neutron shells com- 
pared to the proton shells [sj , and partially to the exis- 
tence of many deformed proton shell closures in the re- 
gion 108 < Z < 120, which drive nuclei with non-magic 
neutron number toward deformation. 
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B. Existence of Shape Isomers 

The existence of a shape (or fission) isomer with defor- 
mation around /? « 1.0 is a prominent feature of actinide 
nuclei [soj . A fission isomer necessarily requires an outer 
fission barrier. The height of the outer barrier with re- 
spect to the isomeric state is shown in Fig. |H1 All forces 
and models confirm the earlier finding that the outer bar- 
rier fades away for transactinide nuclei, see e.g. Ref. [s^ . 
An exception is NL3, which predicts a substantial outer 
barrier for most nuclei, c.f. also figure |H1 As NL-Z2 usu- 
ally does not show a fission isomer, this cannot be a gen- 
eral feature of the RMF model, but has to be a partic- 
ularity of the NL3 parameterization. Similarly, SkI4 is 
an exception among the Skyrme interactions. Most nu- 
clear matter properties of NL3 and NL-Z2 are very close, 
the same holds for the Skyrme interactions SkI4, SkIS 
and SLy6. This suggests that the height of the outer 
barrier is mainly determined by shell structure, not the 
average liquid drop properties. Remember that SkI4 em- 
ploys a non-standard spin-orbit interaction which leads 
to single-particle spectra different from those of the other 
self-consistent models at spherical shape [i^. 



C. Saddle-Point Height 

A most interesting feature of SHE is the stability 
against spontaneous fission. The fission half-live can be 
computed from a tunneling dynamics in the shape de- 
grees of freedom [s^l- Input to that are the collective 
masses along the fission path and the fission barriers. 
We aim here at a mere comparison of stability between 
the different forces. To that end, we assume that the col- 
lective masses are about similar in all cases and confine 
the discussion to the height of the (inner) axial fission 
barrier as key quantity. Keep in mind that the numbers 
given for the barrier height represent an upper limit due 
to possible lowering through triaxial shapes. 

The systematics of the fission barriers for all nuclei 
and mean-field forces considered here are shown in fig- 
ure El We have added results from two other large-scale 
calculations, one employing the semi-classical ETFSI ap- 
proach |30|. and the other within the mic-mac approach 
|2a | . In both cases axial shapes which allow for reflection- 
asymmetry are assumed, similar as in our calculations. 
The mic-mac fission barriers from Ref. |26j are dynami- 
cal barriers, i.e. the barrier which lies on the fission path 
that minimizes the multi-dimensional action. 

All models and forces agree that there is a regime of 
low fission barriers around Z = 110, but the (axial) bar- 
riers increase again when going toward N = 184. There 
are, however, significant differences among the models. 
Comparing the fully self-consistent models, the barriers 
from RMF are much lower than those from SHF. 

And even among the various SHF forces, we see dif- 
ferences in the barriers. It is again SkI4 with its par- 
ticular spin-orbit force which shows the largest barriers. 




164 172 180 164 172 180 
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FIG. 9: The height of the symmetric inner barrier calculated 
in axial symmetry for the models and forces as indicated. We 
have added also ETFSI taken from Ref. 3Q] and mic-mac 
results obtained within the YPE-I-WS model from Ref. j23|. 



Here it is noteworthy that SkI4 predicts largest stability 
for Z = 120 although it places the magic shell closure at 
Z = 114, This shows once again that "magicity" is 
something different from stability as was argued also in 
Ref. 

The ETFSI calculations give barriers that have about 
the size as the ones from the RMF forces, with the differ- 
ence that they produce a bit smaller barriers at the lower 
end and larger regions with more stable nuclei. Similarly, 
the barriers from the mic-mac rather correspond to the 
estimates from the RMF forces. The self-consistent SHF 
results deliver the highest barriers throughout. 

For two nuclei and the forces SLy6 and NL-Z2, the sad- 
dle point shapes are shown in Fig. ^| In all cases, com- 
pact nuclear shapes are obtained. The saddle point defor- 
mation decreases slightly from values around (32 = 0.32 
obtained with both forces for ^^^Hs, to (52 = 0.28 ob- 
tained with SLy6 and (32 — 0.18 predicted by NL-Z2 re- 
spectively for 3°2l20. Note that for NL-Z2 there exists 
no pronounced saddle point, since the barrier is rather 
flat (see Fig.©. 

Figure El also shows fission barriers calculated within 
the mic-mac method |26|. The trends are remarkably 
different from the self-consistent models. A maximum of 
stability is found around Z — 116 and less stability for 
larger systems. Effective mass m* /m = 1 cannot be the 
reason because SkP with m* /m = 1 behaves as all other 
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FIG. 10: Contour plots of the mass density distribution for 
the saddle point configurations of the nuclei ■^*^112i7o (lower 
panels) and ^"^120 (upper panels) as predicted by the forces 
SLy6 (left) and NL-Z2 (right). The z axis is the symmetry 
axis. The contour lines correspond to the densities 0.01, 0.03, 
. . . , 0.15 and 0.17 fm'^ 



SHF and RMF. We can only speculate about possible 
explanations. There might be differences in the smooth 
part of the self-consistent models, e.g. in higher-order 
terms missing in the mic-mac models or even the cur- 
vature term (all of which do not necessarily have the 
right structure in self-consistent models). Or the ef- 
fect comes from the different shell structure between the 
Folded- Yukawa potential used in the mic-mac calcula- 
tions and self-consistent models 01 (although the dif- 
ference among, e.g., SkI4 and NL-Z2 is quite significant, 
while the global trend of the barriers is not). 

It can be speculated that this is a consequence of miss- 
ing shape degrees of freedom, either in missing higher- 
order deformations in the mic-mac method or in the ra- 
dial density distribution. It is well-known that multipole 
moments at least up to ^ = 8 have to be taken into ac- 
count to obtain the full shell effect around 108^8162 
It can be expected that even more shape degrees of free- 
dom are necessary to get the full shell effect for the even 
more complex saddle-point configuration. The standard 
mic-mac models also assume that protons and neutrons 
have the same deformation. An exploration of the con- 
sequences of this constraint in the framework of the self- 
consistent Gogny force is given in Ref. js^. Imposing 
the same deformation for neutrons and protons leads to 
larger barriers of the order of 1 MeV in actinides. A sim- 
ilar effect can be expected for SHE. Even more severe 
might be the parameterization of the radial density dis- 
tribution in mic-mac models. There is no radial degree 
of freedom at all, although it is well-known that chang- 



TABLE II: Comparison of fission barriers (see text). Ex- 
perimental data are taken from Ref. 1S|. The assignment of 
the neutron number has some uncertainty, therefore the same 
experimental barrier appears for two nuclei in the list. 
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ing the surface diffuseness might significantly change the 
ground-state shell correction of SHE. The current pa- 
rameterizations of mic-mac models also prevent "semi- 
bubble" density distributions that might appear at the 
upper end of the nuclei in Fig. El 

Although each missing degree of freedom causes a loss 
in binding energy, the deformation dependence of the var- 
ious effects can be expected to be very different. Depend- 
ing on if the missing energy is larger at the ground state 
or around the saddle point the fission barriers are either 
increased or decreased. This might explain the difference 
in the global trend between mic-mac and self-consistent 
models. 

Table lU compares our calculated barrier heights with 
the lower limits of the barriers of some very heavy nu- 
clei recently deduced from an analysis of the available 
data for fusion and fission 'l3] . Surprisingly these exper- 
imentally estimated barrier heights are similar, or even 
slightly larger, than that of actinide nuclei in the ^"^^Pu 
region. Experimental and calculated values are in agree- 
ment for the Skyrme forces SLy6, SkI3 and SkI4, while 
both RMF forces and the Skyrme interaction SkP signif- 
icantly underestimate the barrier. 

Although the barrier heights are comparable, the life- 
times corresponding to these barriers are much shorter 
than for actinide nuclei as the barriers are much nar- 
rower. For the adjacent 170 110 a fission lifetime of about 
2^1/2 = '''.Blj'g s was reported in Ref. ;8|. This is quite 
short but two orders of magnitude longer than results of 
the mic-mac model which predict about 10~^s |2^. As 
those heavy nuclei are solely stabilized by shell effects 
js^, reliable predictions will be a difficult task for any 
model. 



VII. SEARCH FOR UNDERLYING 
MECHANISMS 

The above results on the fission barriers and its trends 
toward the heaviest SHE show systematic differences be- 
tween SHF and RMF. In this section, we want to ponder 
a bit about possible reasons. For the further discussion it 
is useful to distinguish between the macroscopic part of 



13 




0.0 0.2 0.4 0.6 0.0 0.2 0.4 0.6 

FIG. 11: Inner axial fission barrier (top), proton levels (mid- 
dle) and neutron levels (bottom) for the nucleus ^"^120 with 
NL-Z2 (left) and SLy6 (right). Solid (dotted) lines in the 
Nilsson plots denote single-particle states with positive (neg- 
ative) parity, while the dashed line plots the Fermi energy. In 
the upper panel for NL-Z2 the radial distribution of the total 
density along the z axis is also shown. 

the models (which determines the nuclear matter prop- 
erties and the average trends), and the microscopic part 
(which determines the actual shell structure) . As we will 
see, it is not yet fully clear which part is responsible for 
the observed systematic differences. 



A. Macroscopic aspects 

It is well known that most nuclear matter properties 
from SHF and RMF models differ significantly, see e.g. 
Ref. [^3 and references therein. In lowest order, the bar- 
rier heights can be expected to scale with the surface 
energy coefficient flgurt- With the exception of NL3, how- 
ever, the values for Osurf are quite close for all forces used 
here, see Tabled NL3 gives in most cases smaller (inner) 
barriers in spite of its larger surface energy coefficient. 

If the differences seen in Fig.|Slare rooted in the macro- 
scopic part, another bulk property than Ogurf has to be 



responsible. It is unlikely that this is the volume energy 
because it is basically independent of the nuclear shape. 
The situation is more involved for the volume symme- 
try energy, coefficient Osym. At first glance, the volume 
energy also scales with the nuclear volume, but there is 
an implicit surface effect due to a correlation between the 
neutron skin and asym- the skin increases with increasing 
asym- A systematic study performed in Ref. ^] suggests 
that this relation is unique. A variation of other isovec- 
tor properties as the sum rule enhancement factor ktrk 
or the surface-asymmetry coefficient Kgym leaves the skin 
unchanged. 

And indeed, when we calculate fission barriers with 
the systematically varied Skyrme forces from Ref. [87l| . 
we find that the barriers increase with decreasing Osym, 
while they do not change when the isovector effective 
mass is varied. This is consistent with our findings for 
the barriers, where the Skyrme forces which all have 
Osym ~ 32 MeV have larger barriers than the RMF forces 
with Osym ~ 39 MeV, c.f. tabled This correlation is not 
unique and might apply only to forces fitted accordin g to 
the protocol of Ref. 88]. The Skyrme mass fit MSk7^ 
which follows a very different fitting strategy, has a signif- 
icantly smaller Osym = 27.95 MeV than the Skyrme forces 
used here, but predicts also significantly smaller fission 
barriers for heavy and superheavy nuclei , which are 
in fact similar to our RMF results. With that, a satisfy- 
ing explanation of the difference between RMF and SHF 
concerning fission barrier heights based on macroscopic 
properties of the models is still missing. 



B. Microscopic aspects 

For nuclei with flat or unstable macroscopic poten- 
tial energy surfaces the height of the fission barrier is 
determined by the variation of the shell correction en- 
ergy iSshcii with deformation. i?shcii reflects the deviation 
of the actual density of single-particle levels around the 
Fermi energy ep from an averaged level density. Pivotal 
for the barrier is not the absolute value of i^shoii, but its 
variation, that reflects the change of the single-particle 
spectra with deformation. 

Nilsson plots of the single-particle energies of the nu- 
cleus 3°2l20i82 as calculated with the RMF force NL-Z2 
and the Skyrme interaction SLy6 are shown in figure [TTl 
There are significant differences between the forces which 
can be traced back to their different shell structure at 
spherical shape {^2 = 0) [4^ li^ . The small spin-orbit 
splitting of the 3p and 2/ states obtained with NL-Z2 
leads to a major shell closure at Z — 120, while Z — 120 
is a subshell closure only for SLy6, which (for this neutron 
number) has more prominent gaps in the single-particle 
spectrum at Z = 114 and Z = 126. On the other hand, 
there is a huge gap in the neutron spectrum at A'^ = 184 
for SLy6, while there are several small gaps at = 172, 
182 and 184 for NL-Z2. 

It has to be stressed, however, that for self-consistent 
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models such Nilsson plots cannot be extrapolated very far 
from the N and Z they are calculated for. For deformed 
shapes, the self-consistent optimization of higher multi- 
pole moments when changing N and Z might change the 
single-particle spectra significantly. Additionally, the ra- 
dial shape of the density distribution might change with 
nucleon numbers [s^. 143. l90l| or deformation. Figure ITTI 
provides an example for the latter: at small deforma- 
tion NL-Z2 predicts a semi-bubble shape for ^°^120i82 
(see the small inserts in the upper panel), which around 
(^2 ~ 0.12 changes abruptly into a more regular density 
distribution, thereby causing the discontinuity in the 
single-particle spectra at that deformation. The semi- 
bubble shape reduces the spin-orbit splitting (c.f. the 
proton 2/, 3p or neutron 3d states). The appearance 
of such semi-bubble shapes is force dependent, for SLy6 
it occurs for Z = \2Q only at smaller neutron numbers 
around N — 172. For a more detailed discussion of this 
phenomenon see Refs. 34, 43, 9^. 

Comparing the spectra from NL-Z2 and SLy6 at spher- 
ical shape, there is another significant difference besides 
the spin-orbit splitting: the highly degenerated liii/2+ 
proton and IJ13/2- neutron states above the Fermi en- 
ergy are much lower for NL-Z2 than for SLy6. This does 
not yet lead to a significantly smaller (total) shell cor- 
rection for NL-Z2 (E'shnn = -13.1 MeV) than for SLy6 
[E^^cxi = -14.1 MeV) but their splitting with de- 

formation brings more of these levels close to the Fermi 
energy around the fission barrier for NL-Z2, which might 
be the reason for the difference between NL-Z2 and SLy6. 
From looking at the single-particle spectra in figure 1111 
alone, however, this cannot be decided. To resolve this 
issue, a calculation of the shell correction for deformed 
shapes along the strategy of Ref. [s^ seems highly desir- 
able. 

The single-particle spectra at the spherical point also 
determine the existence or non-existence of the reflection- 
asymmetric outer barrier, and with that of the fission iso- 
mer. At deformations around /32 ~ 0.5, there are several 
single-particle states originating from the intruder and 
the major shells above and below coming close to the 
Fermi energy, see Fig. ^]for the example of '^°^120i82- 
Spectra for other nuclei in this region look quite simi- 
lar. Octupole deformation mixes states with the same 
angular momentum but opposite parity, so depending on 
the actual level ordering it will increase or decrease the 
level density at the Fermi surface. The differences in 
the relative distance of the single-particle states found at 
spherical shape [i^l are reflected in the spectra for these 
very deformed shapes. For example, the subtle shift in 
single-particle energies between the RMF forces NL-Z2 
and NL3 removes the outer asymmetric barrier for the 
former, but not the latter. The single-particle spectra 
at spherical shape including those far above and below 
the Fermi energy, have to be described with very high 
precision to decide if there exist superdeformed states in 
superheavy elements. 




FIG. 12: Single-particle spectra in the superdeformed con- 
figuration obtained from reflection-symmetric calculations for 



VIII. CONCLUSIONS 

We have investigated the systematics of fission bar- 
riers in superheavy elements with Z = 108, . . . , 120 
as predicted by self-consistent mean-field models. As 
typical representatives, we employed the non-relativistic 
Skyrme-Hartree-Fock (SHF) model as well as the rel- 
ativistic mean- field model (RMF), and for each case, 
we used a selection of different parametrizations to ex- 
plore the variances in the predictions. All calculations 
have been done with axial symmetry but allowing for 
reflection-asymmetric shapes. 

As a benchmark for our mean-field models and forces, a 
selection of actinide nuclei ranging from Th to Cf isotopes 
has been utilized to study the predictions for the inner 
and outer axial barriers as well the excitation energies of 
the isomer. Overall, a model dependence of the results 
has surfaced: RMF forces tend to lower, and often too 
low, barriers and excitation energies, while most Skyrme 
forces tend to higher values, which sometimes leads to an 
overestimation on the mean-field level. 

For superheavy nuclei all models and forces agree on 
the systematic trends concerning the fission barriers, 
ground state deformations, and fission modes. There are 
differences in detail. Fission isomers are generally sup- 
pressed in superheavy elements, with the exception of 
the RMF force NL3 in our sample. For larger nuclei, 
there emerge large and systematic discrepancies between 
SHF and RMF concerning the barrier heights, reaching 
a factor two for Z = 120. This amplifies and confirms 
the tendency which has been demonstrated for actinides. 
The factor two in fission barriers around Z — 120 means 
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in absolute numbers that SHF fission barriers are about 
5 MeV larger and this amounts to many orders of mag- 
nitude longer fission lifetimes. Moreover, the SHF forces 
employed here predict also larger barriers than the more 
phenomenological mic-mac models. The reason for the 
systematic difference between SHF and RMF have yet to 
be found out. We suspect that the difference is caused 
by a different shell structure in the models. This point 
deserves more investigation. 

The further systematic trends are shared by SHF 
and RMF. There is a transition from well-deformed 
ground states around Z = 108 to nearly spherical ones at 
Z — 120, which develops through very soft nuclei which 
might exhibit shape isomerism. And there is a marked 
breakdown of fission stability around Z = 110 in agree- 
ment with experimental findings where the a chains of 
superheavy elements are limited at the lower end by fis- 
sion. 

The axially symmetric fission barriers are, of course, 
only a first indicator of fission stability (to be more pre- 
cise an upper limit). One needs yet to include triaxial 
degrees of freedom and to model the dynamics of fission 



to obtain life times. Taking into account the results from 
actinides and superheavy nuclei the mean-fields models 
used in this study seem to deliver lower (and probably 
too low - RMF) and upper (SHF) limits for the barrier 
heights. It is yet difficult to map these differences on spe- 
cial and isolated features of the models which remains an 
urgent and important task for the near future. 
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